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Abstract. Using three parametrizations of the gamma ray burst count data comparison 
is made to cosmological source models. While simple models can fit (V/Vmax) and faint 
end slope constraints, the addition of a logarithmic count range variable describing the 
curvature of the counts shows that models with no evolution or evolution power law in 
redshift with index less than 10 fail to satisfy simultaneously all three descriptors of the 
burst data. The cosmological source density that would be required for a fit is illustrated. 

1. Introduction 

Gamma ray bursts present an astrophysical puzzle as to their mechanism, source, and even 
locality. Despite knowledge of their existence for almost 25 years and over 1500 detected 
bursts, we still do not know if they belong to the realm of the solar system, galaxy, or high 
redshift cosmology, a state that would be unthinkable for other initially mysterious sources 
such as quasars. In the absence of a sufficiently robust clue from their time structure, 
energy spectrum, or a non gamma ray counterpart we turn to their spatial distribution. 
Their angular distribution on the sky is isotropic to 1% and no distance parametrization is 
known for them. We therefore have recourse only to the number-flixx relation, log A^ — log S, 
although for gamma ray bursts the peak count rate in the detector, C, is used instead of 
the energy flux S. 

Section 2 discusses the count distribution from the point of view of properties of 
the detector, including threshold bias, and of the cosmology, e.g. geometry. Section 3 
introduces three parametrizations of the distribution for comparison of theoretical source 
models with the data. The main, new results of this paper are illustrated by Figure 
1 showing the lack of agreement between the data and those cosmological models with 
no evolution or evolution power law in redshift and by Figure 2 showing what form of 
cosmological source distribution would be required to flt the data. 



2. Count distribution 

In analogy with other wavelength regimes (e.g. radio) of astrophysics, one considers the 
data set composed of values of the peak counts C and the count detection threshold 
Ciim- One wants to flnd the number distribution N{C) in the hopes of this giving a clue 
to the source model. (For gamma ray bursts the count rate C is used instead of the 
energy flux 5"). An advantage often accrues to using a relative variable C/Cum, in that 
many instrumental biases cancel out in such a ratio (see, for example, Schmidt, Higdon, 
and Hueter 1988). From this basic relative variable one can form the "volume fraction" 
variable V/Vmax = {C /Ciim)~^^'^ ■ However the effect of variation of the threshold value on 
recovering the entire distribution N{C) is problematic. Another difficulty is the nontrivial 
physical, i.e. spatial, interpretation of variables like this. 
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2.1. Source and detector variation 

First we consider to what extent v = V/Vmax reflects the true source space distribution 
rather than effects arising from variations in source or detector properties. In §2.2 difficul- 
ties from the spatial geometry itself are discussed. 

Define the relation between detected counts C, source count luminosity L, and "dis- 
tance" r according to the usual inverse square law, 

r = (L/47rC)i/2. (1) 

When C is replaced by Cum this likewise defines the maximum "distance" D at which a 
soTircc of (count) luminosity L could be seen by a detector with coimt threshold Cum- Only 
when the source and detector properties, i.e. L and Cum, are known can v — (r/D)^ be 
interpreted in terms of a true distance or spatial measure. A range of L and Cum around 
some fiducial values will cloud our understanding of the true source radial distribution. 
[Note that other variables could have a similar effect on interpreting r as a distance. For 
example, a distribution of source beaming properties, rather than isotropic emission, would 
affect the 47t factor in (1); energy spectrum variations enter the conversion of L from an 
energy luminosity to a count rate; detector efficiency or sky coverage fluctuations alter the 
effective maximum volume, etc.] 

The distribution function p{v) is the fraction of detectable sources lying in a (pseudo) 
volume interval: 

p{v)dv = AfdV/Ntot, (2) 

where A/" is the source number density and Ntot is the total number of detectable sources, 
normalizing the distribution. Given dispersion in a property intrinsic to the source, say L, 
and one intrinsic to the detector, say Cum, how is the distribution function affected relative 
to the standard candle, standard camera case? In other words, how are the distributions 
that include the property dispersions, denoted \p{v;L)]c, [p{v', Ciim)]Lj and \p{v)]l,Cj re- 
lated to the fiducial case p{v; L, CumY- The subscript notation indicates a weighted average 
with respect to that quantity, i.e. an integration over its variation. 

The procedure is most familiar in the case of a luminosity function with J\f = 

n{V)^{L)dL and n the spatial number density. One cannot simply write the integration 
over the dispersion as 

[p(^;CHm)]L = J dL^{L)p{v-Cum.L), (3) 

because of the dependence in the normalizing factor Ntot- There is uniform weighting not 
in terms of the source fraction intrinsically existing at luminosity L, i.e. ^ (-£'), but in the 
fraction detectable., as the definition of Ntot states. That is, 

b(^;C/im)]L = j dL^{L)n{V){dV/dv)l J dL^{L) J dV n{V) 

= J dL^{L)Ntot{L,Cum)p{v,Cum,L) I j dL^L) J dVniV) (4) 

= J dLilj{L,Ciim)p{v;Ciim,L). 
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Two alternative views of this procedure are equivalent: it is either a number weighting [by 
$(L) NtotiL) / Ntot] of the variable p{v; Cum, L) = N^l{L)n{V){dV/dv), as in the last line, 
or a uniform weighting [by ^{L)] of the fundamental variable N^ln{V){dV/dv), as in the 
first line. 

The case of dispersion in the threshold Cum is completely analogous since again this 
variable enters in the normalization factor. Both dispersions together simply involve a 
double integration over the variables. In this case ^{L) is joined by the analogous /{Cum) 
and the weighting function i/j is generalized to 

X{Cum,L) = f{CumML)Ntot{L,Cum)/ J dCumf{Cum) J dL^{L) J dVn{V). (5) 

The effect is to introduce a coarse grained averaging or diffusivity over the distribution 
of V, smearing out sharp features. Even worse, it can act as a redistribution mechanism, 
shifting high count rates (low v), for example, to low. This is very similar to the case of 
amplification bias from gravitational lensing in fiux limited surveys (cf. Schneider 1987). 
The observed faint end slope of —1 in N{C) could be obtained from a steeper true source 
distribution if there existed a broad (fiat /) threshold variation. Petrosian (1993; also see 
Caditz and Petrosian 1993) deconvolves the threshold variations from the count data using 
a nonparametric maximum likelihood statistical method to derive a true faint end slope 
for N{C) in the range —1.6 to —1.8. The distribution is still reasonably fit by a broken 
power law, with a bright end slope of —2.5. See §3 for further discussion. 

2.2. V/Vmax properties 

Because of its historical importance in the gamma ray burst debate we briefiy consider the 
volume fraction variable v = V/Vmax, concentrating on its implications for the geometry 
of the source distribution, particularly cosmologically. The comoving volume element is 

dV ={1 + zfrldrprdio, (6) 

where z is the redshift, rpr the proper distance, the angular diameter distance, and du 
the solid angle element. Since observations indicate that bursts are distributed isotropi- 
cally, duj will simply contribute a factor 47r. 

As known from gravitational lensing, when dealing with a line of sight toward a source 
one must use an angular diameter distance that incorporates the fact that this direction 
is clear of density inhomogeneities that might appreciably (i.e. detectably) lens and alter 
the fiux. These are known as clumpy universe distances and are parametrized by the ratio 
of smoothly distributed matter density to the total (Dyer and Roeder 1973, Ehlers and 
Schneider 1986). When this ratio is unity the angular distances agree with the standard 
Friedmann ones. But calculations indicate that models using smooth vs. totally clumpy 
distances deviate by less than 20% in {V/Vmax), so we will only present results for the 
smooth distances. In general we find that a fiat or open clumpy model acts intermediate 
between the smooth cases f2 = and Q = 1. 
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Given the source functions ^{L) and n{z) and the cosmological model, one can deter- 
mine the distribution and moments of v for comparison with data: 



rV2 i-oo i-z{L,V2) 

/ dvp{v) = NiJ / dL $(L) / dz {drpr/dz)rl{l + zfn{^), 

Jvi Jo Jz{L,vi) 

Ntot= \ dL^iL) \ dzidvprldzylil^zfn^z), (7) 
Jo Jo 

poo rz{L) 

{v) = N,-l / dL $(L) D-\L) / dz {drpr/dz)rlrf{l + zfn{^). 
Jo Jo 

Here z{L, v) is the inversion of ri{z) = D{L)v^^^ with z{L) = z{L, 1) and D following from 
(1) with C = Ciim- The BATSE detector team of the Compton Gamma Ray Observatory 
reports {v) = 0.33 d= 0.01. 

[Note that in contrast to the quasar V/Vmax test (Schmidt 1968), uniformity of source 
spatial distribution does not imply {v) = 1/2 and vice versa. To see this we write the z 
integral of (7) in general functional terms and apply a theorem from integral calculus. 

{v) = 1^ dzj{z)v{z)j ^ dzj{zy, (8) 

in (7) we have j(^z) = {n/Ntot)dV/dz. Now (8) has the form of a weighted average of v{z) 
where v takes the values and 1 at the limits. However, the theorem states that if the 
average of a function with that behavior equals 1/2 over all intervals, i.e. for all values 
Z, then the differential of the function must equal the weighting function, / = dv/dz. 
Using V = (ri/D)^ and equation (6) we see this is not generally true for n uniform, so 
for gamma ray bursts the properties of (f) = 1/2 and uniformity are not equivalent. As 
a specific example, take the model $7 <^ 1, n{z) = no(l + z)"^: one obtains (v) = 1/2 for 
all z, since here coincident ally / = dv/dz. In the quasar test, however, v is defined as 
dV/ /J^'""" dV, so dv/dz = {J^"""^ dV)-'^dV/dz, which does equal / for n uniform.] 
Quasilocal cosmological models, those with detection depth z <ti 1, can be severely 
constrained by calculating the volume relation to first order in redshift: 

(V/V^a.) = (1/2) [1 - (3/7)(l - f3/4)z]. (9) 

Here /3 is the density evolution index in n{z) = no{l + z)^ , i.e. /3 = corresponds to 
no evolution, and this first order result is independent of O. Obtaining (V/Vmax) = 0.33 
would require z — 0.79 or negative evolution, the first violating our low redshift assumption 
and the second requiring excessive evolution at low redshifts (e.g. P = —28 at 2; = 0.1 in 
order to give the value 0.33). 

Were there no problems with threshold bias one could fit simple cosmological models 
with greater depth to the data, obtaining the values (V/Vmax) = (0.4, 0.35, 0.3) at detection 
limits z — (0.7, 1.3, 2.2) for O 1 and ai z — (0.7, 1.5, 3.3) for O 0, with no evolution. 
Excellent approximations to the cosmological model source counts of (7) over the redshift 
range of interest, including evolution, are given to 5% by (v) = 0.5y~^ where y = 1 + z 
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and b = 0.36 - 0.12/3 for = and b = 0.43 - 0.15/3 for ft = 1. Not only the mean but 
the binned distribution (first fine of equation 7) is fit admirably. 



3. Cosmological parametrization 

However, two problems exist with using solely the V/Vmax distribution to describe the 
source counts. One, the information from high counts, including the break away from the 
—5/2 slope, is scrunched into a small region due to transformation to the {C/Ciim)~^^'^ 
variable. Two, the threshold variation problem of Cum also warps the fundamental N(C) 
distribution, as clearly shown by Petrosian's (1993) deconvolution. 
The definition of the differential source counts is 



In the case of uniformity (i.e. sources distributed homogeneously in Euclidean space) the 
volume element dV — Airr'^dr and = r so the integral can be transformed to the form 



Although the bright end slope is matched, there is no possibility of a break in the power law 
with respect to C, regardless of luminosity function. So a luminosity function alone cannot 
give the observed behavior of the bright end slope of —5/2 turning over to a shallower faint 
end slope. There are only four possibilities for the cause of a change in slope. Changes can 
occur by varying the geometry, i.e. dV, by evolution, i.e. n or L{y), and by combinations 
of these with each other or with the luminosity function ^{L). 

Initially neglect the luminosity function; assume all sources have identical luminosity 
L. Including the appropriate delta function for $(!/) in (10) gives 



with the bracketed term evaluated dXy — Y implicitly defined by L = 4:7rCrf{Y). The faint 
slope index becomes constant in the asymptotic, i.e. high redshift, regime with (negative) 
values found to be (6 + /3)/4 or (3 + 2/5)/4 for = or 1. 

The effect of cosmology is that as more distant sources are considered the slope of the 
logA^ — logC relation becomes shallower (greater than —5/2), i.e. there are fewer faint 
sources than expected because of the limited volume available: because of expansion the 
younger universe was smaller. Positive density evolution (more sources at higher redshift: 
P > 0), however, acts to counteract this. At the end of this section we examine how 
unreasonable the density run has to be to match the source counts. Positive luminosity 
evolution increases Y and steepens (shallows) the slope for = (1). The effects are 
opposite because the two models' volume-redshift dependences are very different. 

If the BATSE data are accepted at face value (i.e. without threshold deconvolution), 
the faint end slope is found to be in the range —0.8 to —1.0. In this case low density 




(10) 




(11) 



N{C) = 27rnoC-i [/"Vf (drp^/dr/)] 



(12) 
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cosmological models are poor fits since substantial evolution is needed to achieve such 
a shallow slope (by redshifts of a few rather than asymptotically). Were this extreme 
evolution in fact to exist, however, the source counts would not resemble a broken power 
law but a very gradual turnover. Constraints such as {V /Vmax) from §2.2 add to the 
difficulty, even for = 1 models. 

To put these constraints on a firmer footing we consider three parametrizations of 
the cosmological log — log C distribution for comparison with the data. Two are the 
slope of the differential counts and the logarithmic range of counts between given values 
of the slope. The virtue of using slopes and logarithmic ranges lies in the ability to 
neglect questions of absolute values for the number of bursts or the peak count rate, made 
difficult by detector dependent effects. (Of course once a model achieves a fit under these 
conditions then the normalization is important in recovering the physical values of source 
luminosity and density). We retain {V/Vmax) as a third parameter, both for historical 
reasons and as a complete, rather than point, measure of the distribution. Although these 
are not independent, e.g. knowledge of the slope at every point allows reconstruction of 
the (unnormalized) distribution and hence any of the other measures, there is sufficient 
difference in the areas of emphasis each places on the counts that it is useful to consider 
all three. 

We elaborate the source model to include a version of the cosmological k-correction, 
taking into account the finite energy window through which the bursts are observed and 
the redshift of the source spectrum. The relation between the peak count rate C and the 
received fiux /j^ is 

C = j dvS{v){U/hv), (13) 

where <S'(z^) is the window frequency response and h is Planck's constant. The k-correction 
occurs in transforming the specific luminosity Lj^ to f^i 



L[v{l^z)] d[v{l + z) 



^ = (1 + ^) 4 . - (14) 



where the initial factor of 1 + z accounts for the time dilation since C measures a count 
rate not number. Adopting the spectral models L^, — Lo(i^/t'o)~" and S fiat between 
and 1/+ yields 

C = Lqv^ ^ diyiy-^-''/47Th{iy+ - (1 + zf-''rf'^{z), (15) 

generalizing (1). This differs from the usual inverse square distance law by the 1 + z 
factors; one comes from bandwidth dilation, one from time dilation, and the a dependence 
from the spectral shift through the observing window. The advantage of using differential 
logarithmic parameters is that the entire group of constants in brackets can be neglected. 

Suppose all sources are identical, i.e. same Lq and a, but with evolution possible in 
the source population, as powers of 1 + 2; with /3 the index for density evolution and 7 for 
luminosity {Lq). Generalizing (12), 

N{C) = 47rnoC-'Y^-'+\drpr/dy) \d{y' rf) / dy\-^ = AC-^f{Y), (16) 
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where 5 = 7 + 2 — a and derivatives are evaluated at y = y, the solution to a modified 
equation (1) [or (15)]: C = ^^^-^(y). 

Taking advantage of the form of (16) in solving for the wanted parameters, the (neg- 
ative) slope of the cumulative distribution A^(> C) is 



d\nN{>C) CN{C) 



dlnC N{> C) 



^f{Y)/ J drpry-'rf. (17) 



This can be seen by direct evaluation of J dC C ^f{Y) or by using that A^(> C) is just 
the effective "volume" out to Y{C) to obtain the denominator. The differential slope is 

_ dlnNjC) Cdfdy 

= - dlnC =^-7d^dC- ^^^^ 

Note that because of the bending of the source counts caused by cosmology the relation 
= — 1 docs not generally hold, as it would in the scale free case. 

The volume fraction parameter takes its usual definition in terms of counts relative 
to threshold, {V/Vmax) — {{C/Ciim)~^^'^), but note that it is now even further from 
interpretation as a physical volume than before. Not only is {C/Ciim)~^^'^ not the volume 
element (6) relative to maximum detectable volume (as it would be in the Euclidean case), 
but it is not even the previous distance ratio (ri/D)^ because of the k-correction and time 
dilation factors. 

Finally, define a count range parameter R by 



R{sd) = log[C{sd = 2.25)/C{sd)]. (19) 

This statistic measures how quickly in count space the geometric and evolutionary effects 
become appreciable. It is interpreted as the range to the highest count rate C where N{C) 
has slope Sd from where the slope begins to deviate significantly (10%) from the Euclidean 
value of 5/2. 

By defining the survey depth for each variable - Zy, Zg, and Zr - where the cosmological 
model behavior fits the observations, we discover as viable those models which agree on 
this maxiumum redshift within the data uncertainties. Figure 1 plots the model results 
for f2 = 1 and a variety of evolutionary indices. 

Using data deconvolved of threshold effects by the method of Petrosian (1993) suggests 
parameter values in the ranges of sa = 1.6 — 1.8, R = 1.0 — 1.5, and (independent of 
dcconvolution) (V /Vmax) = 0.33 ± 0.01 . None of the models in Figure 1 simultaneously 
satisfy these values within the errors, as clearly quantified by the depth criterion mentioned 
above. For example, for the (/?, 5) = (—1,0) model to achieve (V /Vmax) = 0.33 ± 0.01 
requires z^ = 1.01 ± 0.1 while Sd = 1.7 ± 0.1 gives Zs = 0.58 ± 0.13 and R = 1.25 ± 0.25 
needs Zr = 0.49 ± 0.13. While the plot is for = 1, lower densities fit even more poorly. 

Recall that the curves incorporate evolution in number density and luminosity, as well 
as k-correction and time dilation, through the parameters P and 6. To obtain a unique 
survey redshift, the curves for (V/Vmax) and Sd need to be shifted to smaller redshifts, 
which can be accomplished by decreasing P or 5; each affects the results in roughly the 
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same magnitude. A unique depth, within the data uncertainties, is not achieved until 
f3 + S < —10, at which point the depth has decreased to 2; ~ 0.1. That requires an extreme 
amount of evolution in a very small redshift range. 

The count range parameter proves very useful, giving a measure of the rate of cur- 
vature of the count distribution. Without paying heed to this constraint simple, even 
nonevolutionary, models will seem to fit for depths varying from z ^ 1 — 7. However the 
logarithmic count ranges are then some two times larger than the values indicated by the 
deconvolved data, i.e. a factor of 60-250 off the observed count rate. 

Of course one can always adjust the source spatial distribution in (10) so as to obtain 
any count behavior desired. Figure 2 illustrates the spatial density n{z) of gamma ray 
burst sources required to match the observed count distribution. We adopt two forms for 
this distribution: a broken power law with bright slope -2.5 and faint slope -1.7, and a more 
gradual behavior that interpolates between these two asymptotes. For the cosmology we 
take models with = and Q — 1 and incorporate luminosity evolution, k-correction, and 
time dilation effects through the parameter 5 as before. The requisite cosmological density 
behavior for explanation of the observations is not satisfied by the form n{z) ~ (1 + -2)^, 
which would be a straight line on the plot with slope /3, or by any familiar cosmological 
source distribution. 

4. Conclusion 

The gamma ray burst source count plot is similar to the classical Hubble diagram in that 
(for the cosmological scenario) the curvature away from the Euclidean behavior is a direct 
measure of cosmological effects. The bright end slope of —5/2 has a natural interpretation 
in terms of nearby sources, z < 0.1, where geometry and evolution have not yet broken 
source uniformity. The faint slope value determines a characteristic survey depth, or 
redshift, as does the weighted average (V/Vmax) of the entire data set, and the observed 
range in counts between the transition region and the asymptotic behavior. 

These three distinct estimates for the depth out to which gamma ray bursts are 
detected - Zgiope, Zy, and Zrange ^ must be equal for the cosmological model fit. Figure 1 
shows, however, that this does not hold for any of the models considered with power law 
evolution. Figure 2 illustrates the run of spatial number density required for cosmological 
sources to match the observed count distribution, again in disagreement with power law 
evolution. The models took into account the cosmological characteristics of geometry, 
k-correction and time dilation, as well as number and luminosity evolution power law in 
redshift. 

Acknowledgments. I am indebted to Dieter Hartmann, Wlodek Kluzniak, and Vahe Pet- 
rosian for helpful discussions. Use was made of the BATSE online database at gronews® 
gr ossc . gsfc . nasa. gov. 
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Figure captions 

Fig. 1. The statistical parameters, normahzed to the center of their observed range, are 
plotted vs. redshift. Those three curves starting a,t z — with value 1.52 show (V/Vmax), 
those with value 1.47 show s^, and the positive slope curves show R. Solid curves have 
5)=(0,0), long dashed have (0,-1), and short dashed (-1,0). The intersection of each 
curve with the value one defines the appropriate depth z^, Zg, or z^. A successful model 
would have all solid curves, say, intersecting the dotted line simultaneously. Of course 
uncertainty in the data broadens the dotted line, into 1 ± 0.03 for (V/Vmax), 1 =t 0.06 for 
Sd, and 1 ± 0.2 for R. 

Fig. 2. The cosmological spatial number density derived to match a source counts model 
is plotted vs. redshift. The solid and dotted lines use a broken N{C) and $7 = and 1, 
respectively; the short and long dashed lines a gradual N{C) and = and 1, respectively. 
The break, or turnover, in counts is assumed to lie at 2 = 0.1; it is actually a function of 
the intrinsic luminosity and spectrum and will scale the horizontal axis accordingly, but 
should lie close to the value adopted. 
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